New tissue priors for improved automated classification of subcortical brain structures on MRI.

Despite the constant improvement of algorithms for automated brain tissue classification, the accurate delineation of subcortical structures using magnetic resonance images (MRI) data remains challenging. The main difficulties arise from the low gray-white matter contrast of iron rich areas in T1-weighted (T1w) MRI data and from the lack of adequate priors for basal ganglia and thalamus. The most recent attempts to obtain such priors were based on cohorts with limited size that included subjects in a narrow age range, failing to account for age-related gray-white matter contrast changes. Aiming to improve the anatomical plausibility of automated brain tissue classification from T1w data, we have created new tissue probability maps for subcortical gray matter regions. Supported by atlas-derived spatial information, raters manually labeled subcortical structures in a cohort of healthy subjects using magnetization transfer saturation and R2* MRI maps, which feature optimal gray-white matter contrast in these areas. After assessment of inter-rater variability, the new tissue priors were tested on T1w data within the framework of voxel-based morphometry. The automated detection of gray matter in subcortical areas with our new probability maps was more anatomically plausible compared to the one derived with currently available priors. We provide evidence that the improved delineation compensates age-related bias in the segmentation of iron rich subcortical regions. The new tissue priors, allowing robust detection of basal ganglia and thalamus, have the potential to enhance the sensitivity of voxel-based morphometry in both healthy and diseased brains.

Age- and sex-related variations in the brain white matter fractal dimension throughout adulthood: an MRI study.

PURPOSE: To observe age- and sex-related differences in the complexity of the global and hemispheric white matter (WM) throughout adulthood by means of fractal dimension (FD). METHODS: A box-counting algorithm was used to extract FD from the WM magnetic resonance images of 209 healthy adults from three structural layers, including general (gFD), skeleton (sFD), and boundaries (bFD). Model selection algorithms and statistical analyses, respectively, were used to examine the patterns and significance of the changes. RESULTS: gFD and sFD showed inverse U-shape patterns with aging, with a slighter slope of increase from young to mid-age and a steeper decrease to the old. bFD was less affected by age. Sex differences were evident, specifically in gFD and sFD, with men showing higher FDs. Age × sex interaction was significant mainly in the hemispheric analysis, with men undergoing sharper age-related changes. After adjusting for the volume effect, age-related results remained approximately the same, but sex differences changed in most of the features, with women indicating higher values, specifically in the left hemisphere and boundaries. Right hemisphere was still more complex in men. CONCLUSIONS: This study is the first that investigates the WM FD spanning adulthood, treating age both as a continuous and categorical variable. We found positive correlations between FD and volume, and our results show similarities with those investigating small-world properties of the brain networks, as well as those of functional complexity and WM integrity. These suggest that FD could yield a highly compact description of the structural changes and also might inform us about functional and cognitive variations.

Disentangling in vivo the effects of iron content and atrophy on the ageing human brain.

Evidence from magnetic resonance imaging (MRI) studies shows that healthy aging is associated with profound changes in cortical and subcortical brain structures. The reliable delineation of cortex and basal ganglia using automated computational anatomy methods based on T1-weighted images remains challenging, which results in controversies in the literature. In this study we use quantitative MRI (qMRI) to gain an insight into the microstructural mechanisms underlying tissue ageing and look for potential interactions between ageing and brain tissue properties to assess their impact on automated tissue classification. To this end we acquired maps of longitudinal relaxation rate R1, effective transverse relaxation rate R2* and magnetization transfer - MT, from healthy subjects (n=96, aged 21-88 years) using a well-established multi-parameter mapping qMRI protocol. Within the framework of voxel-based quantification we find higher grey matter volume in basal ganglia, cerebellar dentate and prefrontal cortex when tissue classification is based on MT maps compared with T1 maps. These discrepancies between grey matter volume estimates can be attributed to R2* - a surrogate marker of iron concentration, and further modulation by an interaction between R2* and age, both in cortical and subcortical areas. We interpret our findings as direct evidence for the impact of ageing-related brain tissue property changes on automated tissue classification of brain structures using SPM12. Computational anatomy studies of ageing and neurodegeneration should acknowledge these effects, particularly when inferring about underlying pathophysiology from regional cortex and basal ganglia volume changes.

Multivariate decoding of brain images using ordinal regression.

Neuroimaging data are increasingly being used to predict potential outcomes or groupings, such as clinical severity, drug dose response, and transitional illness states. In these examples, the variable (target) we want to predict is ordinal in nature. Conventional classification schemes assume that the targets are nominal and hence ignore their ranked nature, whereas parametric and/or non-parametric regression models enforce a metric notion of distance between classes. Here, we propose a novel, alternative multivariate approach that overcomes these limitations - whole brain probabilistic ordinal regression using a Gaussian process framework. We applied this technique to two data sets of pharmacological neuroimaging data from healthy volunteers. The first study was designed to investigate the effect of ketamine on brain activity and its subsequent modulation with two compounds - lamotrigine and risperidone. The second study investigates the effect of scopolamine on cerebral blood flow and its modulation using donepezil. We compared ordinal regression to multi-class classification schemes and metric regression. Considering the modulation of ketamine with lamotrigine, we found that ordinal regression significantly outperformed multi-class classification and metric regression in terms of accuracy and mean absolute error. However, for risperidone ordinal regression significantly outperformed metric regression but performed similarly to multi-class classification both in terms of accuracy and mean absolute error. For the scopolamine data set, ordinal regression was found to outperform both multi-class and metric regression techniques considering the regional cerebral blood flow in the anterior cingulate cortex. Ordinal regression was thus the only method that performed well in all cases. Our results indicate the potential of an ordinal regression approach for neuroimaging data while providing a fully probabilistic framework with elegant approaches for model selection.

PRoNTo: pattern recognition for neuroimaging toolbox.

In the past years, mass univariate statistical analyses of neuroimaging data have been complemented by the use of multivariate pattern analyses, especially based on machine learning models. While these allow an increased sensitivity for the detection of spatially distributed effects compared to univariate techniques, they lack an established and accessible software framework. The goal of this work was to build a toolbox comprising all the necessary functionalities for multivariate analyses of neuroimaging data, based on machine learning models. The "Pattern Recognition for Neuroimaging Toolbox" (PRoNTo) is open-source, cross-platform, MATLAB-based and SPM compatible, therefore being suitable for both cognitive and clinical neuroscience research. In addition, it is designed to facilitate novel contributions from developers, aiming to improve the interaction between the neuroimaging and machine learning communities. Here, we introduce PRoNTo by presenting examples of possible research questions that can be addressed with the machine learning framework implemented in PRoNTo, and cannot be easily investigated with mass univariate statistical analysis.

Recommendations to improve imaging and analysis of brain lesion load and atrophy in longitudinal studies of multiple sclerosis.

Focal lesions and brain atrophy are the most extensively studied aspects of multiple sclerosis (MS), but the image acquisition and analysis techniques used can be further improved, especially those for studying within-patient changes of lesion load and atrophy longitudinally. Improved accuracy and sensitivity will reduce the numbers of patients required to detect a given treatment effect in a trial, and ultimately, will allow reliable characterization of individual patients for personalized treatment. Based on open issues in the field of MS research, and the current state of the art in magnetic resonance image analysis methods for assessing brain lesion load and atrophy, this paper makes recommendations to improve these measures for longitudinal studies of MS. Briefly, they are (1) images should be acquired using 3D pulse sequences, with near-isotropic spatial resolution and multiple image contrasts to allow more comprehensive analyses of lesion load and atrophy, across timepoints. Image artifacts need special attention given their effects on image analysis results. (2) Automated image segmentation methods integrating the assessment of lesion load and atrophy are desirable. (3) A standard dataset with benchmark results should be set up to facilitate development, calibration, and objective evaluation of image analysis methods for MS.

Regional specificity of MRI contrast parameter changes in normal ageing revealed by voxel-based quantification (VBQ).

Normal ageing is associated with characteristic changes in brain microstructure. Although in vivo neuroimaging captures spatial and temporal patterns of age-related changes of anatomy at the macroscopic scale, our knowledge of the underlying (patho)physiological processes at cellular and molecular levels is still limited. The aim of this study is to explore brain tissue properties in normal ageing using quantitative magnetic resonance imaging (MRI) alongside conventional morphological assessment. Using a whole-brain approach in a cohort of 26 adults, aged 18-85years, we performed voxel-based morphometric (VBM) analysis and voxel-based quantification (VBQ) of diffusion tensor, magnetization transfer (MT), R1, and R2* relaxation parameters. We found age-related reductions in cortical and subcortical grey matter volume paralleled by changes in fractional anisotropy (FA), mean diffusivity (MD), MT and R2*. The latter were regionally specific depending on their differential sensitivity to microscopic tissue properties. VBQ of white matter revealed distinct anatomical patterns of age-related change in microstructure. Widespread and profound reduction in MT contrasted with local FA decreases paralleled by MD increases. R1 reductions and R2* increases were observed to a smaller extent in overlapping occipito-parietal white matter regions. We interpret our findings, based on current biophysical models, as a fingerprint of age-dependent brain atrophy and underlying microstructural changes in myelin, iron deposits and water. The VBQ approach we present allows for systematic unbiased exploration of the interaction between imaging parameters and extends current methods for detection of neurodegenerative processes in the brain. The demonstrated parameter-specific distribution patterns offer insights into age-related brain structure changes in vivo and provide essential baseline data for studying disease against a background of healthy ageing.

Do clients with acquired brain injury use the splints prescribed by occupational therapists? A descriptive study.

Clients with acquired brain injury often demonstrate hypertonicity and decreased function in their upper limbs, requiring appropriate intervention. Splinting is one of the intervention methods that is widely used to address these issues. Literature shows that some clients are not using splints following fabrication. However, there is a paucity of research about the factors that influence clients to use or not use splints. This study aims to investigate these influential factors for clients with upper limb hypertonicity. Two survey tools including therapist and client questionnaires were developed and completed by both therapists and clients. Six therapists and 14 clients participated in this study and completed the relevant questionnaires. The results illustrate that most clients (13 out of 14) were continuing to use their splints four weeks following discharge from hospital. The main goals of choosing splints for both therapists and clients were prevention of contracture and deformity. The most indicated client reasons for adhering to the splint wearing program were therapist-related factors including clients' trust and reliance on their therapists. Further reasons for clients implementing the recommended splint-wearing program and clinical implications are discussed.

Genotype-phenotype interactions in primary dystonias revealed by differential changes in brain structure.

Our understanding of how genotype determines phenotype in primary dystonia is limited. Familial young-onset primary dystonia is commonly due to the DYT1 gene mutation. A critical question, given the 30% penetrance of clinical symptoms in DYT1 mutation carriers, is why the same genotype leads to differential clinical expression and whether non-DYT1 adult-onset primary dystonia, with and without family history share pathophysiological mechanisms with DYT1 dystonia. This study examines the relationship between dystonic phenotype and the DYT1 gene mutation by monitoring whole-brain structure using voxel-based morphometry. We acquired magnetic resonance imaging data of symptomatic and asymptomatic DYT1 mutation carriers, of non-DYT1 primary dystonia patients, with and without family history and control subjects with normal DYT1 alleles. By crossing the factors genotype and phenotype we demonstrate a significant interaction in terms of brain anatomy confined to the basal ganglia bilaterally. The explanation for this effect differs according to both gene and dystonia status: non-DYT1 adult-onset dystonia patients and asymptomatic DYT1 carriers have significantly larger basal ganglia compared to healthy subjects and symptomatic DYT1 mutation carriers. There is a significant negative correlation between severity of dystonia and basal ganglia size in DYT1 mutation carriers. We propose that differential pathophysiological and compensatory mechanisms lead to brain structure changes in non-DYT1 primary adult-onset dystonias and DYT1 gene carriers. Given the range of age of onset, there may be differential genetic modulation of brain development that in turn determines clinical expression. Alternatively, a DYT1 gene dependent primary defect of motor circuit development may lead to stress-induced remodelling of the basal ganglia and hence dystonia.

Automatic detection of preclinical neurodegeneration: presymptomatic Huntington disease.

BACKGROUND: Treatment of neurodegenerative diseases is likely to be most beneficial in the very early, possibly preclinical stages of degeneration. We explored the usefulness of fully automatic structural MRI classification methods for detecting subtle degenerative change. The availability of a definitive genetic test for Huntington disease (HD) provides an excellent metric for judging the performance of such methods in gene mutation carriers who are free of symptoms. METHODS: Using the gray matter segment of MRI scans, this study explored the usefulness of a multivariate support vector machine to automatically identify presymptomatic HD gene mutation carriers (PSCs) in the absence of any a priori information. A multicenter data set of 96 PSCs and 95 age- and sex-matched controls was studied. The PSC group was subclassified into three groups based on time from predicted clinical onset, an estimate that is a function of DNA mutation size and age. RESULTS: Subjects with at least a 33% chance of developing unequivocal signs of HD in 5 years were correctly assigned to the PSC group 69% of the time. Accuracy improved to 83% when regions affected by the disease were selected a priori for analysis. Performance was at chance when the probability of developing symptoms in 5 years was less than 10%. CONCLUSIONS: Presymptomatic Huntington disease gene mutation carriers close to estimated diagnostic onset were successfully separated from controls on the basis of single anatomic scans, without additional a priori information. Prior information is required to allow separation when degenerative changes are either subtle or variable.

Atrophy progression in semantic dementia with asymmetric temporal involvement: a tensor-based morphometry study.

We performed a longitudinal anatomical study to map the progression of gray matter atrophy in anatomically defined predominantly left (LTLV) and right (RTLV) temporal lobe variants of semantic dementia (SD). T1-weighted MRI scans were obtained at presentation and one-year follow-up from 13 LTLV, 6 RTLV, and 25 control subjects. Tensor-based morphometry (TBM) in SPM2 was applied to derive a voxel-wise estimation of regional tissue loss over time from the deformation field required to warp the follow-up scan to the presentation scan in each subject. When compared to controls, both LTLV and RTLV showed significant progression of gray matter atrophy not only within the temporal lobe most affected at presentation, but also in the controlateral temporal regions (p<0.05 FWE corrected). In LTLV, significant progression of volume loss also involved the ventromedial frontal and the left anterior insular regions. These results identified the anatomic substrates of the previously reported clinical evolution of LTLV and RTLV into a unique 'merged' clinical syndrome characterized by semantic and behavioral deficits and bilateral temporal atrophy.

Voxel-wise analysis of diffusion tensor MRI improves the confidence of diagnosis of corticobasal degeneration non-invasively.

Corticobasal degeneration (CBD) presents with symptoms that often overlap with other neurological conditions. In many cases, diagnosis, prognosis and consequent clinical management remain uncertain. Structural and functional asymmetric brain changes represent the most consistent imaging findings that may assist in CBD diagnosis. Diffusion Tensor MRI (DT-MRI) is a quantitative technique that allows microscopic tissue abnormalities to be non-invasively assessed in vivo. A single case of clinically suspected CBD with symmetric diffuse brain atrophy on conventional-MRI scans was studied using DT-MRI by voxel-wise comparison with eight healthy subjects. The lateralized distribution of DT-MRI abnormalities was consistent with clinical features providing a substantial support to the diagnosis.

Diffusion-based spatial priors for imaging.

We describe a Bayesian scheme to analyze images, which uses spatial priors encoded by a diffusion kernel, based on a weighted graph Laplacian. This provides a general framework to formulate a spatial model, whose parameters can be optimized. The application we have in mind is a spatiotemporal model for imaging data. We illustrate the method on a random effects analysis of fMRI contrast images from multiple subjects; this simplifies exposition of the model and enables a clear description of its salient features. Typically, imaging data are smoothed using a fixed Gaussian kernel as a pre-processing step before applying a mass-univariate statistical model (e.g., a general linear model) to provide images of parameter estimates. An alternative is to include smoothness in a multivariate statistical model (Penny, W.D., Trujillo-Barreto, N.J., Friston, K.J., 2005. Bayesian fMRI time series analysis with spatial priors. Neuroimage 24, 350-362). The advantage of the latter is that each parameter field is smoothed automatically, according to a measure of uncertainty, given the data. In this work, we investigate the use of diffusion kernels to encode spatial correlations among parameter estimates. Nonlinear diffusion has a long history in image processing; in particular, flows that depend on local image geometry (Romeny, B.M.T., 1994. Geometry-driven Diffusion in Computer Vision. Kluwer Academic Publishers) can be used as adaptive filters. This can furnish a non-stationary smoothing process that preserves features, which would otherwise be lost with a fixed Gaussian kernel. We describe a Bayesian framework that incorporates non-stationary, adaptive smoothing into a generative model to extract spatial features in parameter estimates. Critically, this means adaptive smoothing becomes an integral part of estimation and inference. We illustrate the method using synthetic and real fMRI data.

The role of the medial temporal lobe in autistic spectrum disorders.

The neural basis of autistic spectrum disorders (ASDs) is poorly understood. Studies of mnemonic function in ASD suggest a profile of impaired episodic memory with relative preservation of semantic memory (at least in high-functioning individuals). Such a pattern is consistent with developmental hippocampal abnormality. However, imaging evidence for abnormality of the hippocampal formation in ASD is inconsistent. These inconsistencies led us to examine the memory profile of children with ASD and the relationship to structural abnormalities. A cohort of high-functioning individuals with ASD and matched controls completed a comprehensive neuropsychological memory battery and underwent magnetic resonance imaging for the purpose of voxel-based morphometric analyses. Correlations between cognitive/behavioural test scores and quantified results of brain scans were also carried out to further examine the role of the medial temporal lobe in ASD. A selective deficit in episodic memory with relative preservation of semantic memory was found. Voxel-based morphometry revealed bilateral abnormalities in several areas implicated in ASD including the hippocampal formation. A significant correlation was found between parental ratings reflecting autistic symptomatology and the measure of grey matter density in the junction area involving the amygdala, hippocampus and entorhinal cortex. The data reveal a pattern of impaired and relatively preserved mnemonic function that is consistent with a hippocampal abnormality of developmental origin. The structural imaging data highlight abnormalities in several brain regions previously implicated in ASD, including the medial temporal lobes.

Progression of structural neuropathology in preclinical Huntington's disease: a tensor based morphometry study.

BACKGROUND AND OBJECTIVES: Regional cerebral atrophy occurs in carriers of the Huntington's disease (HD) gene mutation before clinical diagnosis is possible. The current inability to reliably measure progression of pathology in this preclinical phase impedes development of therapies to delay clinical onset. We hypothesised that longitudinal statistical imaging would detect progression of structural pathology in preclinical carriers of the HD gene mutation, in the absence of measurable clinical change. METHODS: Thirty subjects (17 preclinical mutation positive, 13 mutation negative) underwent serial clinical and magnetic resonance imaging (MRI) assessments over an interval of 2 years. Statistically significant changes in regional grey and white matter volume on MRI were analysed using tensor based morphometry (TBM). This technique derives a voxel-wise estimation of regional tissue volume change from the deformation field required to warp a subject's early to late T1 images. RESULTS: Over 2 years, there was progressive regional grey matter atrophy in mutation-positive relative to negative subjects, without significant clinical progression of disease. Significant grey matter volume loss was limited to bilateral putamen and globus pallidus externa (GPe), left caudate nucleus, and left ventral midbrain in the region of the substantia nigra. CONCLUSIONS: While these results are consistent with previous cross sectional pathologic and morphometric studies, significant progression of atrophy in HD before the onset of significant clinical decline is now demonstrable with longitudinal statistical imaging. Such measures could be used to assess the efficacy of potential disease modifying drugs in slowing the progression of pathology before confirmed clinical onset of HD.

The precision of anatomical normalization in the medial temporal lobe using spatial basis functions.

We investigated the accuracy of spatial basis function normalization using anatomical landmarks to determine how precisely homologous regions are colocalized. We examined precision in terms of: (1) the number of nonlinear basis functions used by the normalization procedure; (2) the degree of (Bayesian) regularization; and (3) the effect of substituting different templates and how this interacted with the number of basis functions. The face validity of spatial normalization was assessed as a function of these parameters, using the colocalization of homologous landmarks in a test sample of 20 normally developing children and 5 children with bilateral hippocampal pathology. Our results suggest that when optimal normalization parameters are used, anatomical landmarks in the medial temporal lobes are colocalized to within a standard deviation of about 1 mm. When suboptimal parameters are used this standard deviation can increase up to 3 mm. Interestingly the optimal parameters are those that provide a rather constrained normalization as opposed to those that optimize intensity matching at the expense of rendering the warps "unlikely." The implications of our results, for users of voxel-based morphometry, are discussed.

Distributional assumptions in voxel-based morphometry.

In this paper we address the assumptions about the distribution of errors made by voxel-based morphometry. Voxel-based morphometry (VBM) uses the general linear model to construct parametric statistical tests. In order for these statistics to be valid, a small number of assumptions must hold. A key assumption is that the model's error terms are normally distributed. This is usually ensured through the Central Limit Theorem by smoothing the data. However, there is increasing interest in using minimal smoothing (in order to sensitize the analysis to regional differences at a small spatial scale). The validity of such analyses is investigated. In brief, our results indicate that nonnormality in the error terms can be an issue in VBM. However, in balanced designs, provided the data are smoothed with a 4-mm FWHM kernel, nonnormality is sufficiently attenuated to render the tests valid. Unbalanced designs appear to be less robust to violations of normality: a significant number of false positives arise at a smoothing of 4 and 8 mm when comparing a single subject to a group. This is despite the fact that conventional group comparisons appear to be robust, remaining valid even with no smoothing. The implications of the results for researchers using voxel-based morphometry are discussed.

Classical and Bayesian inference in neuroimaging: theory.

This paper reviews hierarchical observation models, used in functional neuroimaging, in a Bayesian light. It emphasizes the common ground shared by classical and Bayesian methods to show that conventional analyses of neuroimaging data can be usefully extended within an empirical Bayesian framework. In particular we formulate the procedures used in conventional data analysis in terms of hierarchical linear models and establish a connection between classical inference and parametric empirical Bayes (PEB) through covariance component estimation. This estimation is based on an expectation maximization or EM algorithm. The key point is that hierarchical models not only provide for appropriate inference at the highest level but that one can revisit lower levels suitably equipped to make Bayesian inferences. Bayesian inferences eschew many of the difficulties encountered with classical inference and characterize brain responses in a way that is more directly predicated on what one is interested in. The motivation for Bayesian approaches is reviewed and the theoretical background is presented in a way that relates to conventional methods, in particular restricted maximum likelihood (ReML). This paper is a technical and theoretical prelude to subsequent papers that deal with applications of the theory to a range of important issues in neuroimaging. These issues include; (i) Estimating nonsphericity or variance components in fMRI time-series that can arise from serial correlations within subject, or are induced by multisubject (i.e., hierarchical) studies. (ii) Spatiotemporal Bayesian models for imaging data, in which voxels-specific effects are constrained by responses in other voxels. (iii) Bayesian estimation of nonlinear models of hemodynamic responses and (iv) principled ways of mixing structural and functional priors in EEG source reconstruction. Although diverse, all these estimation problems are accommodated by the PEB framework described in this paper.

Classical and Bayesian inference in neuroimaging: applications.

In Friston et al. ((2002) Neuroimage 16: 465-483) we introduced empirical Bayes as a potentially useful way to estimate and make inferences about effects in hierarchical models. In this paper we present a series of models that exemplify the diversity of problems that can be addressed within this framework. In hierarchical linear observation models, both classical and empirical Bayesian approaches can be framed in terms of covariance component estimation (e.g., variance partitioning). To illustrate the use of the expectation-maximization (EM) algorithm in covariance component estimation we focus first on two important problems in fMRI: nonsphericity induced by (i) serial or temporal correlations among errors and (ii) variance components caused by the hierarchical nature of multisubject studies. In hierarchical observation models, variance components at higher levels can be used as constraints on the parameter estimates of lower levels. This enables the use of parametric empirical Bayesian (PEB) estimators, as distinct from classical maximum likelihood (ML) estimates. We develop this distinction to address: (i) The difference between response estimates based on ML and the conditional means from a Bayesian approach and the implications for estimates of intersubject variability. (ii) The relationship between fixed- and random-effect analyses. (iii) The specificity and sensitivity of Bayesian inference and, finally, (iv) the relative importance of the number of scans and subjects. The forgoing is concerned with within- and between-subject variability in multisubject hierarchical fMRI studies. In the second half of this paper we turn to Bayesian inference at the first (within-voxel) level, using PET data to show how priors can be derived from the (between-voxel) distribution of activations over the brain. This application uses exactly the same ideas and formalism but, in this instance, the second level is provided by observations over voxels as opposed to subjects. The ensuing posterior probability maps (PPMs) have enhanced anatomical precision and greater face validity, in relation to underlying anatomy. Furthermore, in comparison to conventional SPMs they are not confounded by the multiple comparison problem that, in a classical context, dictates high thresholds and low sensitivity. We conclude with some general comments on Bayesian approaches to image analysis and on some unresolved issues.